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We give an extensive treatment of the pairing symmetry in the ferromagnetic superconductor 
UGei. We show that one can draw important conclusions concerning the superconducting state, 
considering only the transformation properties of the pairing function, without assumptions about 
the form of the pairing amplitudes. 



I. INTRODUCTION 

The interplay of superconductivity and magnetism is a 
subject of central importance to accomplish a basic un- 
derstanding of the fundamental interactions in a solid. In 
particular, the coexistence of superconducting and ferro- 
magnetic states is still a matter of controversy. Early 
contributions to this topic include the works of Matthias 
and Bozorth on ZrZn^ and of GinzburgS, who was the 
first to address the problem of coexistence of both phe- 
nomena. 

It was only recently that UGe^r (for a more complete 
reference see^), was found to display such properties. 
Since magnetic fields tend to suppress conventional su- 
perconductivity, such a coexistence is expected to take 
place in a unconventional stated. 

Although one can find many efforts in the search for 
a microscopic theory specific to this problem (— ~— ), there 
is still a lot of work to be done, and a phenomenological 
approach is of great importance. 

The theory for unconventional superconductivity was 
given ir£ for systems where the normal state (non super- 
conducting) is paramagnetic, and its group of symmetry 
includes the time reversal operation. In turn, for ferro- 
magnetic metals, the normal state is associated to a mag- 
netic group, or cogroup, and the superconducting states 
arising from this normal state are classified according to 
the co- representations of this cogroup 9 . 

Consequences for ZrZri2 were explored, including 
both, strong and weak spin orbit couplings^. General 
discussions on other systems were also given ir>ii. Fol- 
lowing the above approaches, in this paper we propose a 
detailed and rigorous treatment of this problem for the 
case of UGe2, in the strong spin orbit regime. 



II. THE BASIC HAMILTONIAN AND 
SYMMETRIES. 



Our Hamiltonian is based on some aspects of the 
physics of UGe^r^- Its ferromagnetic (FM) properties 
resemble closely a metallic state, where the Fermi surface 
has a big splitting (70mev) between the up and down spin 
bands. This fact points for a choice of an 'equal spin pair- 
ing' (ESP) state, with intraband pairing only. Moreover, 
the FM state is strongly anisotropic. The magnetization 
is clearly pinned along the crystallographic a axis, which 
we take as the z direction. This is very important, since 



the transverse magnetic fluctuations are pair breaking 
for our proposed ESP state. Hence, these fluctuations 
are considered to be much smaller than the interactions 
along the z axis, and we keep just this component. 

The basis of our theory is a generalized BCS Hamil- 
tonian that allows for triplet paring—. We also add a 
term describing a ferromagnetic exchange coupling with 
strength J— . After a standard mean field approximation, 
we remain with: 
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with the definitions of the mean fields M and A aa (fc) 
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M acts as a potential describing the magnetization of 
the system and A CTCr (fc) is the pair potential of the super- 
conducting state. In our case, this potential is indeed a 
matrix that we write in the form: 

A _ ( A(fc) tt \ 

which can be parametrized in terms of the elegant and 
useful d(k) vector formalism^: 



i(d(k).a)a y = 



-d x (fc) + id y (k) d z {k) 



d z (k) 



d x (fc) 4- id y (k) 



(4) 

With the choice of an ESP state, we are left with a two- 
dimensional order parameter lying on the xy plane, while 
the spin of the Cooper pairs is along the z axis: 



d x {k) 
d v {k) 



(A(k) u - A(fc) tt ) 
(A(fc) n + A(fc) u ). 
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Table I: Corepresentations for the Cogroup D 2 (C 2 ). 



r 


E 


C 2z 


KC 2x 


KC 2 y 


.4 


1 


1 


1 


1 


B 


1 


-1 


1 


-1 



For our purposes, a more suitable d(k) vector is: 

d r (k) = ±[-(x + iy)A(k) n + (x-iy)A(k) u ] (5) 

where the T index labels the irreducible corepresentation 
associated with the transition. This d(k) vector is the 
order parameter of our theory for triplet pairing. We 
denote by G the normal (ferromagnetic) state from which 
the superconducting state will rise. It is written as: 

G = M x 17(1) , (6) 

where M is the magnetic group of the system, con- 
structed from its crystal symmetry, and U(X) is the gauge 
group. UGe2 has an orthorhombic crystal structure with 
full inversion symmetry^, identified with the C mmm space 
group. From now on, we treat only the point group 
symmetry, as it is common when one consider the Lan- 
dau theory of phase transitions. It is based on the pre- 
scription that macroscopic properties in thermodynamic 
equilibrium should not depend on details related to the 
translational symmetry of the system. The point group 
of interest is D 2 h, which is generated by the operations 
{C 2z , C 2x , I}. To describe a state with a magnetization 
vector along the z axis, one needs to include the time- 
inversion operator K, combined with some symmetry op- 
eration. A suitable choice that solves this problem is: 

{E,C 2zi KC 2x ,I} (7) 

yielding the complete magnetic point group: 

M = {E, C 2z , KC 2x , KC 2y , 1} . (8) 

Note that a magnetization (axial vector) along the z axis 
is completely invariant under the operations of as it 
should be. To further simplify our work, we only pro- 
ceed with the group {E,C 2z ,KC 2x ,KC 2x }, taking care 
in fixing the parity of the representations. This cogroup 
is denoted D^iC^), and has two irreducible corepresen- 
tations, which we label A and B (see table U). 

From table [I] we see that Sa is associated with the 
trivial breaking of U(l) gauge symmetry, while for Sb 
we should consider a subgroup of ?7(1), isomorphic to 
the invariant subgroup {E,KC 2x \. Such a subgroup is 
simply {E, exp(i7r)}. The result is that Sb contains the 
so-called non-trivial elements. Table [TTI summarizes these 
results. 



III. ANALYSIS OF RESULTS 

Each of these superconducting classes is related to a 
gap structure. This is determined from the transforma- 



Table II: Allowed Superconducting Classes. 



Classes 


Elements 


Sa(D 2 (C 2 )) 


{E,C 2z ,KC 2x ,KC 2y } 


S B ({E,KC 2x }) 


{E,C 2z e^,KC 2x ,KC 2y e^} 



tion properties of the d(k) vector under the action of the 
class elements. For the C 2z element of the class Sa, and 
taking the k vector as k = (0, 0, k), we have: 

C 2 J A {k ) = -d A {C^k ) = -d A {k Q ), 

but from table U we get: 

C 2 J A (k ) = di A (k ) =>• ^(fco) - , 

meaning that the gap will have point nodes on the inter- 
sections of the Fermi surface with the k z axis. For the Sb 
class, taking C 2z e m and the k- vector as fco = (k x ,k y ,0), 
we have 

C 2z e i ^d B (k ) = ACA) = A-fco) • 

Recalling that A(fc) must satisfy the Pauli antisymmetry 
property A(— k) = — A(fc), and looking at tabled we 
write: 

d*(*0) = C 2z e"d B (k ) = (-l)d B (k ) , 

resulting that the gap will have a line node on the inter- 
section of the Fermi surface with the fc^fc^-plane. Finite 
values for the interband paring potential change this sce- 
nario. To investigate its consequences, we directly study 
the transformation properties of the paring potentials. 
First, we set up the transformation for the Cfe,ct opera- 
tors under (j8]): 

(C 22 |0) : A4 )± -> TiA4, fc)± 

(KC 2x \6) : A4 >± -> ±iA*c! Ca , fci± , (9) 

where A is any complex number. To deduce these rules, 
one should represent the time inversion as K = (io~ v )9 
(where 8 is the complex conjugate operation), and con- 
sider that pseudo-spin states transform like spin eigen- 
states. The translational part of the operation is explic- 
itly written, because we are now dealing with properties 
of microscopic quantities in a way that lattice transla- 
tions are important. The situation is not more involved, 
since the space group of interest is symorphyc. Now, we 
require invariance of the Hamiltonian ([1]), including all 
the A Q ^ potentials in the two body interaction part. Be- 
low, we work out an example for the non-unitary operator 

{Kc 2x \oy. 

^[(KC 2:c |0)A tt (fc)c[. t ][(KC 2 ,|0) C t _ fct ] = 

k 

= £A^(fc)[H) clc^HH) 4^ fct ] 

k 
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= ^A| t (C 2 - 1 A : ')[4, t ][c t _ fc , t ] 
k' 

(KC 2x \0)A n (k) = A* n (C^k) 

We turn now to the expansion of the paring amplitudes in 
terms of the cobasis functions Ar(fc) = XL Vr,ifi(k). The 
?7r,i are the order parameter components usually associ- 
ated with the expansion of the Ginzburg-Landau free en- 
ergy, and fi(k) are the orbital basis functions. We choose 
the phase as KC2 X fr(k) = fr(k) for the basis functions, 
being this choice dictated by the character table of this 
group. This choice implies that we have KCi x f]T = Vr- 
We write 

A aa . B ( k ) = i> B f A {k) , (10) 
with a = J [, 4, and 

A tUA (k) = M A (k) 

A n , B (k) = yj B f B (k) . (11) 

Now, we can see that the symmetry dictated nodes for 
the interband and intraband paring amplitudes are given 
by different corepresentations . If one investigates, for 
example, the gap structure in the (k x ,k y ) plane, even if 
the superconducting state happens to be in the B state 
(realized trough A A ), the above predicted nodal line will 
be absent, since A^-j-j, relates to the f A {k) amplitudes, 
which do not present this line. But, as already discussed, 
we expect that a large splitting of the spin bands hinders 
such interband components. Another interesting feature 
is that the appearance of order parameters of the same 
symmetry in different sheets of the Fermi surface, is con- 
nected with terms like c^^_ k ^.c\, , cL k i^ which do not con- 
serve spin. This possibility is allowed due to the strong 
spin orbit effect present in those systems. 



To investigate further the gap structure, we consider 
specific basis functions which come from the lattice 
Fourier series, thus taking into account the translational 
orthorhombic symmetry^. Considering odd functions, 
we have: 

f(k)=J2 c n sin (k -Rn) (12) 



where R n are the nearest neighbors sites of the UGe2 
Bravais lattice^. Using table 2 we get: 

f(k) A = as in(^)cos(^) (13) 
f(k) B = / 3sin(^)cos(^) (14) 

Clearly, the basis (|13til4[) yields nodes other than the 
ones predicted before ('accidental nodes'). Note that all 
the symmetry dictated nodes are here, and we have just 
determined extra nodes due to this possible particular 
choice of the basis functions. Experiments should be 
used in order to distinguish among the possible forms 
for the paring amplitudes. Theoretical predictions ( 10 i 13 ) 
may allow a meaningful interpretation of experimental 
results. From those, ultrasound attenuation experiments 
are most relevant for this kind of system. Such an ex- 
periment, acting as a microscopic probe, would be very 
important in order to select possible scenarios for the 
f(k) functions. Considerations regarding the structure 
dictated by a particular choice of basis functions are im- 
portant also, as it may act as a guide for the construction 
of a microscopic theory. In principle, thermodynamic and 
transport properties may distinguish between the A and 
B phases, since the topology of the nodes implies differ- 
ent low temperature behaviors for such quantities 5 . 

This work was supported by CNPQ and FAPESP 
(Brasil). 



1 B. T. Matthias & R. M. Bozorth. Phys Rev. 109, 604 
(1957). 

2 V. L. Ginzburg. Soviet Physics JETP 4, 2 (1957). 

3 S. S. Saxena et al. Nature 406, 10 (2000). 

4 A. Huxley et al. Phys Rev B 63, 144519 (2001). 

5 M. Sigrist & K. Ueda. Rev. Mod. Phys. 63, 2 (1991). G. 
E. Volovik and L. P. Gorkov. Sov Phys JETP 61, 4 (1985). 
E. I. Blount. Phys. Rev. B 32, 5 (1985). 

6 T. R. Kirkpatrick, D. Belitz. Phys. Rev. B 67, 024515 
(2003). 

7 Y. Zhou, C. D. Gong. Europhys Lett, 74, 145-150 (2006). 



8 N. Karchev. Phys Rev B 67, 054416 (2003). 

9 V. P. Mineev. Phys Rev B 66, 134504 (2002). 

10 K. V. Samokhin and M. B. Walker. Phys Rev B 66, 024512 
(2002); ibid 174501 (2002). 

11 A. Huxley et al. Physica C 403, 9-14 (2004). 

12 J. B. Ketterson & S.N. Song. Superconductivity. Cambridge 
University Press (1999). 

13 M.B. Walker, M.F. Smith and K. V. Samokhin, Phys. Rev. 
B 65, 014517 (2001). 



